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1. Introduction 

In the present paper, we consider integrable discretizations of the nonhnear partial differential 
equation 

WTXX - 2k^wx + 2wxwxx + wwxxx =0, (1) 
which belongs to the Harry-Dym hierarchy H] 121 [51. Here K is a real parameter and, as shown 
subsequently, can be normalized by the scaling transformation when K 7^ 0. A connection 
between Eq.([T]) and the sinh-Gordon equation was established in |4|. When K = 0, Eq.([Tli is 
called the Hunter-Saxton equation and is derived as a model for weakly nonlinear orientation 
waves in massive nematic liquid crystals [5J. The Lax pair and bi-Hamiltonian structure 
were discussed by Hunter and Zheng 161. The dissipative and dispersive weak solutions were 
discussed in details by the same authors fT^JSl. 

Equation ([U can be viewed as a short-wave model of the Camassa-Holm equation ||9l 

wt + 2k^wx — wtxx + 3wwx — Iwxwxx + wwxxx ■ (2) 
Following the procedure in ||T0lim[T2l . we introduce the time and space variables T and X 

T = eT, X=e-^X, 

where e is a small parameter Then w is expanded as w = £^(wo + £wi + •••) with w,- 
(; = 0, 1 , • • •) being functions of f and X. At the lowest order in e, we obtain 

^o.Txx - 2k^vi'o,x + ^Wq x^o xx + ^^qWo xxx = , (3) 
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which is exactly Eq.([T]) after writing back into the original variables. Based on this fact, 
Matsuno obtained the A^-cuspon solution of Eq.([T]l by taking the short-wave limit on the A^- 
soliton solution of the Camassa-Holm equation [TS^ T?!. 

Note that the parameter K of Eq.([T]) can be normahzed to 1 under the transformation 

which leads to 

Wtxx — 2w.r + 2WxWxx + WWxxx = . (4) 

We call Eq.(|4| the short wave model of the Camassa-Holm equation (SCHE). Without loss 
of generality, we will focus on Eq. (|4]l and its integrable discretizations, since the solution of 
Eq.([T]i with arbitrary nonzero K, its integrable discretizations and the corresponding solutions 
can be recovered through the above transformation. 

The reminder of the present paper is organized as follows. In section 2, we reveal a 
connection between the SCHE and the bilinear form two-dimensional Toda-lattice (2DTL) 
equations. The parametric form of A^-cuspon solution expressed by the Casorti determinant 
is given, which is consistent with the solution given in |13|. Based on this fact, we propose 
an integrable semi-discrete analogue of the SCHE in section 3, and further its integrable full- 
discrete analogue in section 4. The concluding remark is given in section 5. 



2. The connection with 2DTL equations, and A^-cuspon solution in determinant form 

2.1. The link of the SCHE with the two-reduction of2DTL equations 

In this section, we will show that the SCHE can be derived from the bilinear form of two- 
dimensional Toda lattice (2DTL) equations 



(5) 



where is the Hirota D-derivative defined as 

and D_i and D\ represent the Hirota D derivatives with respect to variables x_i and x\, 
respectively. 

It is shown that the A^-soliton solution of the 2DTL equations (|5]l can be expressed as the 
Casorati determinant lfT6l[T7l 
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automatically satisfies the above dispersion relations. 

Applying the two-reduction T„_i = {J^=ipjY^'^n+i^ i-e., enforcing pi — — ^„ / — 
1, • • • ,A^, we get 

)t„-t„=t2^i, (8) 



-D,Dy-\]g-g = f. (10) 



where the gauge transformation i„ — > iJ^^iPi)"ln is used. Letting 1q— /, Ti —g and x_ i —s, 
XI — y, the above bilinear equation (O takes the following form: 

-(^-D,Dy-\^f-f = g\ (9) 
2 

Introducing u ~ gj f, Eqs.(|9]l and ( fTOl i can be converted into 

-(ln/),. + l-M2^ (11) 
-(lng),,+ l (12) 

Subtracting Eq.(fT2b from Eq.lfTTTi. one obtains 

|(lnp)3, + l=p2, (13) 

by letting p = m^. 

Introducing the dependent variable transformation 

w = -2(lng),,, 

it then follows 



(lnp), = -^w,, (14) 



or 

P 

2 

by differentiating Eq.(fT2b with respect to i 
In view of Eq.(fT4li. Eq.(fT3]) becomes 

2 V2 

Introducing the hodograph transformation 



-^V^^y] +1=P • (15) 



x^2y-2{\a.g)s, 
t = s. 



and referring to Eq.(fT2]l. we have 



which implies 
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Thus, Eqs.(fT4li and ( fTSl ) can be cast into 
C id, + wd^)\np 



y -vv„ + 1 = p- 

By eliminating p, we arrive at 



(16) 



{dt + wdx) In (- Wxt + 1 ) = -2wx , 



or 



{dt + wdx)wxx - 2wx (1 - Wxx) = , 

which is actually the SCHE 



2.2. The N-cuspon solution of the SCHE 

Based on the link of the SCHE with the two-reduction of 2DTL equations, the A^-cuspon 
solution of the SCHE (HJi is given as follows: 

w = -2(lng),,, 

x = 2y-2(ln^),, 
t^s, 

,<■•<.' 

l<ij<Af 

V(/p) = + ai2{-piye-P''''-"'y+'^o, . (17) 

Moreover, the A^-cuspon solution of the SCHE ([T]i with non-zero K is given as follows: 
w(y,r) = -2(lng),,, (18) 



(19) 



where 



with 



\<i,i<N 



W = ai^yp'leP'y+'IP'+'^-^+ai2{-pi)"e-''-y-''IP'+<'^ . 



We remark here that to assure the regularity of the solution, the T-function is required to be 
positive definite. In what follows, we list the one-cuspon and two-cuspon solutions. For 
N = I, the T-function is 

by choosing a\,\l a\2 — —\, which yields the one-cuspon solution 

w{y, T) = - ^sech2 [p,{y + kT / p\ +ya)\ , 
Pi 

^ - 1^ - ^ { 1 + tanh [pi (y + Kr/p2 +3^0)] } . 
The profiles of one-cuspon with K = 1 .0 and K = 0. 1 are plotted in Fig. [T] 
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Figure 1. Plots for one-cuspon solution for pi = \/2 and different K: (a) K = 1 .0; (b) K = 0. 1 . 



The T-function corresponding to the two-cuspon solution is 



P\~P2 



with 



9; = 2p, {y + kT/p^ + y/o) , / = 1 , 2 . 
Here ai^i/a\;i — — I and a2, 1/^2,2 = 1 are chosen to assure the regularity of the solution. 



3. Integrable semi-discretization of the SCHE 

Based on the link of the SCHE with the two-reduction of 2DTL equations clarified in the 
previous section, we attempt to construct the integrable semi-discrete analogue of the SCHE. 
Consider a Casorati determinant 
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with \|/) satisfies the following dispersion relations 
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where Ak is defined as Ak\\f{k) = 'fW-j(^-i) ^ jjj particular, we can choose as 

V\"\k) = pni-ap.)->^e^-+q'l{l-aqi)->^e^- , 

. I . 1 
Q = —s + q,o , r|; = —s + y]io , 
Pi <li 



(20) 
(21) 
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which automatically satisfies the dispersion relations (|20] i and ( I^TT i. The above Casorati 
determinant satisfies the bilinear form of the semi-discrete 2DTL equation (the Backlund 
transformation of the bilinear equation of the 2DTL equation) ifTTlfTSl 

^D, - 1 ^ T„ (fc + 1 ) • T„ (A:) + T„+ 1 (A: + 1 )t„_ 1 (fc) = . (22) 

Applying a two-reduction condition pi — —qi, i — 1, • • • ,A^, which implies T„_i ^ in+i, 
we obtain 

-A -l)fk+i- fk = gk+\gk , (23) 
a J 

1 



fl V fk 
1 A 8k+i 



Ds-l j gk+i ■ gk = fk+\gk , (24) 

by letting To(^) ^ fk, li{k) = gk- 

Letting Uk ~ gk/fk, Eqs.(l23]l and (l24l l are equivalent to 

+1 = (25) 

-- (In^ ) +1 =M^-^Mi'. (26) 

Subtracting Eq.(l26b from Eq.(l25Tl. one obtains 

Uk+lUk f, Uk+i\ 2 2 ^^^^ 

——[^^— I +^^''k+l'^k- (27) 
a \ / ^ 

Introducing the discrete analogue of hodograph transformation 



and 



It then follows from Eq. 



Xk = 2ka-2{\ngk)s: 



6/, = xj^^ 1 — x^. = 2fl — 2 ( In ^^ii 

gk 



~ 2a 
Ok = 



Uk+lUk 



or 



4a^ 
h 

by assuming pi. = u^. 

Introducing the dependent variable transformation 

Wk = -2{lngk)ss, 

Eq.jZTli becomes 

- In^ +1_ 0. (29) 

5^ V Pk J s H 



Differentiating Eq.(l26]l with respect to s, we have 

^{wk+\ - Wk) = ^- — {\nuk+iuk)^ = — ^- {\npk+\9k), 

la Uk+iUk lUk+\Uk 
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or 



(lnPA+lPir)s = -^{wk+l~Wk)- 

OA- 
Eliminating Pit and ^k+\ from Eqs.(|29]l and ( |30] l, we obtain 



-^{wk+i-wk) - ^ [wk-Wk-i] = -{hk + hk-\)~'^a — 



or 



5/1-1 



1 



1 



5/t 5i:_l 



7 - / c- 4fl^ 

A = —M Ok — ^ 



(30) 



(31) 



(32) 



by defining a difference operator A and an average operator M as follows 



Fk+\ 



Fk+i +Fk 



(33) 



Furthermore, a substitution of Eq. (l28] l into Eq. (|30] | leads to 
d8k 

-r- = Wk+i-Wk- 
ds 

Equations dSTT i and ( |33] ) constitute the semi-discrete analogue of the SCHE. 

Next, let us show that in the continuous limit, a {5k ^ 0), the proposed semi-discrete 
SCHE recovers the continuous SCHE. To this end, Eqs.dsTli and (l33l l are rewritten as 



■{Awk-Awk-i) + l 



, ds5k = wk+i -wk- 
By taking logarithmic derivative of the first equation, we get 
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The dependent variable w is regarded as a function of x and f , where x is the space coordinate 
of the A;-th lattice point and t is the time, defined by 



k-l 
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t^s. 



In the continuous limit, a — (5it 0), we have 
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ds = 9, + Tr-3,v ^ 3, + wd^ , 
as 

where the origin of space coordinate xq is taken so that ^r-^ cancels wq. Thus the above 

as 

semi-discrete SCHE converges to 

{dt+wdx){-Wxx+'^) 



-Wxx + 1 



or 



(34) 



(3, + wdx)wxx = 2w.t {-Wxx + 1) , 
which is nothing but the SCHE ©. 

In summary, the semi-discrete analogue of the SCHE and its determinant solution are 
given as follows: 

The semi-discrete analogue of the SCHE 

1 2/^1 
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The determinant solution of the semi-discrete SCHE 

Wk = -2{lngk)ss, 

s o fk+lfk 

Ok = Xk+i -Xk = 2a 

gk+igk 

Xk = 2ka - 2{lngk)s , 
t — s , 
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gk 
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, fk 
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l<i.j<N ' 



^Vy\k)=a,,lp■j{l~ap.)-'eP^' ■'+'^°-+a,2{-piy{l+api)-'e-P-' '+"^01 .(36) 

Introducing new independent variables Xk ~ Xk/K and T ~ t/x, we can include the 
parameter K in the semi-discrete SCHE d35T l 



^{Wk+I-Wk) - -^^(wk-Wk-l) = 7r9(§i+S't-l) 

Ok Ok-l 2k^ 



2a' { ^ + ^ 

.Ok Ok-l 



(37) 



Wk+l - Wk , 



d8k 
{ 'df 

where hk = Xk+i — Xk and s = kT. This is the semi-discrete analogue of the SCHE ([T]l. 

The A^-cuspon solution of the semi-discrete SCHE ( [37T i with the parameter K is given by 

Wk = -2{lngk)ss, 

s V V ^'^ fk+lfk 
Ok — ^k+l —^k — , 

K gk+lgk 

Xk='-^-l{\ngk)s, 



l<i.j<N 



7 fk 



VW(^) = fl,up/(l - ap.)-'eP--'-^+'^»-+a,2i-p.yii + apt)-'e-P^''^+< 



l<i.j<N 



(38) 
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4. FuU-discretization of the SCHE 

In much the same way of finding the semi-discrete analogue of the SCHE, we seek for its 
full-discrete analogue and in the process we arrive at its A^-cuspon solution. 
Consider the following Casorati determinant 



l<i,j<N 



(39) 



A-1 ]8k+i,r8kj+fk+ijfkj=0, (43) 



where 

^\f\"\k,l) = ai^ip'lil-api)-'' {l~bpry'e^-+aip_q'l{l^aq,)-'' {l-bqry'e'^- , 
with 

= pr^s+t,io, r[i = qr^s+i\iQ. 

It is known that the above determinant satisfies bilinear equations ifTsl 

QD,-l^T„(^+l,/)-T„(fe,/)+T„+i(fe+l,/K-i(^,0=O, (40) 

and 

{bD, - l)x„{k,l + 1) •T„+i(^,/) +T„(fe,/)T„+i(^,/ + 1) = 0. (41) 

Here a,b aie mesh sizes for space and time variables, respectively. 

Applying the two-reduction T„_i = {Ylf=i pj)^^'^n+i, i-e., enforcing pi = —q,, i = 
1,---,N, and letting i(){k,l) ~ fkj, = gkj, the above bilinear equations take the 

following form: 

-D, - 1 ) fk+i,rfk,i +8k+i,igk,i =0. (42) 
a I 

1 

— J 

a 

{bD, - i -gki+fkjgk.i+i =0, (44) 

(JbD, - l)gu+x -fu +gk,ifk.i+i =0, (45) 

where the gauge transformation T„ (HjLi Pi)"%i is used. It is readily shown that the above 
equations are equivalent to 

« V fu /.I fk+ljfkj 

}-L^J±ll] ^i^fjllilhl^ (47) 
a V gk,l J , gk+l.lgkj 

bLl^\ ^i^fMiMli^ (48) 

V 8kJ J, fkj+lgk,l 

t(lni^) =i^Skjfk,i+i_ ^49^ 

V /*:,/ Js gk.l+\fk,l 

We introduce a dependent variable transformation 

w,,, = -2(lng,,,)^^, (50) 
and a discrete hodograph transformation 

Xkj = 2ka - 2(lngkj)s , (51) 
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then the mesh 

hkj =xk+\i-xk,i^2a-2{\n^^^^\ (52) 



gkj 



is naturally defined. It then follows 



\n^J±ll\ =2«-i(8,,, + 5,_i,,). (53) 



In view of Eq.(|47]i. one obtains 

fk+ljfk.l _ ^kj_ 

8k+i.i8k.i 2a 
A substitution into Eq.(|46]l yields 



fk.l J s 



it then follows 



Starting from an alternative form of Eq.(l47Ti 



we obtain 



'k,i 2fl-2 In^^ V 8k+\J8k,i 



^k-l.l V 8kJ8k-l.l 

Subtracting Eq.(|59b from Eq.dSSTl. we obtain 



Si./ h-l.l \ fk-lj J s \ 8k-lJ 

By using the relations (l53T l and (|56] |. we finally arrive at 



(8^,/ + 8^-i,/)+2fl^ — + - =0. (61) 



(54) 



In — — ) fl - — , (55) 

'k,l 



In^) ^2a-2fl2f 1 + 1 ). (56) 

Jk-ljJ, \OkJ Ok- 1,1^ 



2a-2(ln^) ^2a^^±i^, (57) 

8k,l Js 8k+\,l8k,l 



-2 (In 



Wk+i,t-Wkj \ fikj ),, fk+i.ifkj I ^^g^ 



by taking logarithmic derivative with respect to s. A shift from ^ to ^ — 1 gives 

Wk,l-Wk-l,l f fk.lfk-l,l\ 

In . (59) 



Wk+l,l-WkJ WkJ-Wk-l.l f^^ fk+lj \ A ^^Q^ 



5a-,/ 5/.-1,/ 2 \Okj Ok-ij, 

Similar to Eq.(|32]i. Eq.(|6T]i constitutes the first equation of the full-discretization of the SCHE, 
which can be cast into a simpler form: 

AV,/ = ^Mf5,,/-^y (62) 

Next, we seek for the second equation of the fuU-discretization. Recalling (|46]|-(|49]|. one 
could obtain 

2fl-2fln^^±H±i') fln^^) -i 

Xk+lJ+l —Xk,l+l _ V Js _ V fk+hl J s ^ ^g^-j 

Xk+u~xk,i 2fl-2fln^') fln^) -i 
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here a shift from / to / + 1 in WT\ and a shift from ^ to + 1 in ( |49] l are employed. 
From Eqs.jSOb, ( fSSl l and ( |58] l. one can find the following two relations 

, gk+\,i+\\ Wk+u-Wk,i-2a^ 1 

In— 7 = T5 l-7(-«<:+i,;+^jt,z-2.xi:+i /+ij , (64) 

Ik+\,i J s ^^k,i 4 

, fkj+\\ Wk+\j+\-wu+i+2a^ 1 

In = ^ -:{xk+i.i+i+xk^i+\~2xkj) , (65) 

gk.i /j 25(t^/+i 4 

after some tedious algebraic manipulations. Substituting these two relations into ( |63] ), we 
finally obtain the second equation of the fully discrete analogue of the SCHE 

§*:,/+l -§*:,/ 1 s / , ox 
-, r -:0k J+\ 1 ~2Xki 

b 4 

= ^ (w;t+ 1 ,/+ 1 + 1 - Wk,i+\ -Wk,i) ■ (66) 
Taking the continuous limit — > in time, we have 

b ds ' 

5<r,/+l {xk+\J+\ ~2xi.j) — > 0, 

5/t,/+ 1 5i:,/ 1 ^/ + XkJ - 2XA:+ 1 ^/+ 1 ) , 

and ^ 

2 (v^it+iJ+i + Wk+i,t - Wkj+i - Wk,i) Wk+i - Wk- 

Therefore, one recovers exactly the second equation of the semi-discrete SCHE ([33l l. 

In summary, the fully discrete analogue of the SCHE and its determinant solution are 
given as follows: 

The fully discrete analogue of the SCHE 

s s ::{0k,i + 0k-i,i)+2a [^ + ^ ^0, 

Ok.i Ok-i.i 2 \bk,i Ok-ij J 

^^^^ — ^ + {xk+ij+i +Xk,i+i ~ 2xk,i) (67) 

b 4 

+ ^hAxk+ij+Xk.i-2xk+ij+i) = ^{wk+ij+i +Wk+ij - Wkj+i - Wkj) ■ 
The determinant solution of the fully discrete SCHE 

X ^ h.i8k.i-hl, 

Wkj = -2{\ngk,ij,-. = -2 J , 

8k.i 

Xk I = 2ka — 2{\ngk /).? = 2ka — 2—^ , 

gk,i 

s T fk+l,lfk,l 

Ok.i =Xk+ij—Xkj = la , 

gk+\.lgk.l 
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gk,i 



hi = 



ds 



l<i,j<N 
(0) 



7 fk,l — 



l<i,j<N 



^vf\k,i} ^^ff\k,l) ^^ff\k,l) 
^vf\k,l) ^vf\k,l) 



,.(2) 



.(3), 



hu = 



32 



gk.i 



ds^ 



^v[''\k,l) ^vf\k,l) ^vf\k,l) 



.(3), 



^v[-'\k,l) ^^f^^\k,l) 



.(3), 



• ¥f (^,0 



V)^'^(^,0 ¥}v'(^:0 



„(2) 



.(3), 



••• ¥^(^,0 
^i^fiKl) 4'\k,l) xv'i\k,l) ... 






,.(1) 



.(3), 



y\,l'\k,l)=atjpiil-api)-''{l-bpi-^) ' e^^ + ai^i-piY (1 + apt)-'' {l + bpr^) 'e^' , 

= pr^s + ^,0 , ri; = -pr^s + r|,o . (68) 

Note that s is an auxiliary parameter. By virtue of s, hk.i and Jikj can be expressed 
hk,t — dsgkj and hki — d-^gkj, respectively, because the auxiliary parameter s works 

elements of the above determinant by 3f¥|"'' {k, ^ ¥;" ^\k, 0- 

Introducing new independent variables X^ i = x^ i/k and b = b/K, we can include the 
parameter K in the full-discrete SCHE ( |67] l: 

Wk+lJ~WkA WkJ-Wk-lJ 



as 
on 



5a.,/ 
5it,;+i — Sir,/ 



Ok-l.l 



i 9 / i 

^ (5a:,/ + 5a:_ 1 ,/) + 2fl ( — 



^k.l 



>k-l,l 



' ^;j^5i(:,/+i {Xk+ij+i +Xkj+i — 2Xkj) 



(69) 



+ ^-j^5i:,/ (Xir+i,/ +Xi:,/ —2X^+1, /+l) = - {Wk+lj+l +Wk+lJ —WkJ+l — Wkj) ■ 

Similarly, the A^-cuspon solution of the full-discrete SCHE ( |69l l with the parameter K is 
given as follows: 



Wkj = -2{ln gkj),, = -2 



hk.igkj - hjj 



gii 



2ka 2 2ka 2 h^ / 
Xkj = (ln^/t,/)i = - 

K ' K KgkJ 

s V V 2a fk+\,ifk.i 

Ok.l = Xk+ ij—Xkj — , 

^ gk+i.igk.i 



gkJ = 



l<i,j<N 



fkJ = 



l<i,j<N 
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ds 



ds^ 



¥^''(^,/) 



Mff\k,l) Mf^\k,l) 



^vf\k,l) ¥^(-^,0 



W'}^\k,l) 



xvl^\k,l) xv^^\k,l) 



¥f (fc,/) 



-/ 



= P,' + ^lO , II; = -pi + r|,o . 



(70) 



5. Concluding remarks 



In the present paper, bilinear equations and the determinant solution of the SCHE are obtained 
from the two-reduction of 2DTL equations. Based on this fact, integrable semi-and full- 
discrete analogues of the SCHE are constructed. The A^-soliton solutions of both continuous 
and discrete SCHEs are formulated in the form of the Casorati determinant. Note that the 
short pulse equation was also obtained from the two-reduction of the 2DTL equation 119|. 

Finally, we remark that the present paper is one of our series of work in an attempt 
of obtaining integrable discrete analogues for a class of integrable nonlienar PDEs whose 
solutions possess singularities such as peakon, cuspon or loop soliton solutions. New discrete 
integrable systems obtained in this paper, along with the semi-discrete analogue for the 
Camassa-Holm equation fT5\ and the semi-discrete and fully discrete analogues of the short 
pulse equation 1 19 1 deserves further study in the future. 
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